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FOREWORD 


A better understanding of the basic mechanisms of sediment 
transport by wave action is being sought by the Coastal Engineering 
Research Center (formerly the Beach Erosion Board) through several 
different approaches. One of these which has been pursued under 
contract with the University of California at Berkeley, is reported 
on herein. This report discusses the mechanisms of sediment trans- 
port in a layer immediately adjacent to the ocean floor for long 
waves of small amplitude and relatively deep water. Proceeding on 
the premise that a bed-load function exists for the assumed type 
of flow, a bed-load transport equation was developed. Effects of 
both the unsteady mean flow velocity and the turbulent fluctuations 
are taken into account. From this equation the rate per unit width 
at which sediment in the bed layer is shifted by the oscillatory 
flow is calculated. 


This report was prepared at the Hydraulic Engineering Labora- 
tory of the University of California at Berkeley in pursuance of 
contracts DA-49-055-eng-17 and DA-49-055-CE-63-4 with the Beach 
Erosion Board (now the Coastal Engineering Research Center). These 
contracts provided in part for the study of sand movement by wave 
action. The author of the report, George Kalkanis, was at that 
time a graduate student and research assistant at the University. 


This report is published under authority of Public Law 166, 
79th Congress, approved July 31, 1945, as modified by Public Law 
88-172, approved November 7, 1963. 
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TRANSPORTATION OF BED MATERIAL DUE TO WAVE ACTION 
by 


George Kalkanis, University of California 


ABSTRACT 


A practical method has been developed which can be used to 
determine the rate of sediment transportation in a layer adjacent 
to the ocean floor. The method is applicable only when the flow 
in this layer caused by surface waves is unstable. The waves in 
question should be of small amplitude and great length, permitting 
the linearization of the equations of motion. The fundamental 
principle of the supporting theory is that at equilibrium the sub- 
merged weight of the solid particle resting on the ocean floor is 
balanced by the vertical component of the resultant of all the 
hydrodynamic forces exerted on the particle by the flow above. 
Both effects of the unsteady mean flow velocity as well as of the 
turbulent fluctuations are taken into account. The distribution 
of the lift forces associated with the former was determined 
experimentally while a statistical approach based on the experi- 
ence with the same phase of the problem in a steady mean flow was 
used to determine the latter. Proceeding on the premise that a 
bed-load function exists for the type of flow we are dealing with 
here the bed-load equation was developed. From this equation it 
is possible to determine the rate at which sediment in the bed- 
layer is shifted by the oscillatory flow across a section of unit 
width, The concentration of bed material in the layer that at any 
instance is at a state of oscillatory motion can be determined 
from the above rate. Finally, this concentration in combination 
with the velocity distribution in the bed-layer associated with 
any incidental secondary flow is used to calculate the rate of 
transportation of bed material in the direction of this flow. 


1, INTRODUCTION 


The purpose of this study is to develop a method by means of which 
low rates of sediment transport due to the action of surface waves may be 
predicted with sufficient accuracy. More specifically the waves considered 
are long waves of small amplitudes in relatively deep water. The problem 
is of great importance to engineers and scientists as can be witnessed 
from the considerable amount of research devoted to it especially since 
the beginning of the century. An extensive survey of the literature, how- 
ever, revealed that so far the subject has been treated only qualitatively. 


The numerous publications are mostly restricted in presenting observa- 
tions made in the field or the laboratory. With exception of a few 
purely empirical equations, which describe only particular measurements, 
no attempt has been made, at least to the writer's knowledge, to derive 
basic relationships which will constitute the basis of a quantitative 
treatment, In this paper we will try to analyze separately the indivi- 
dual mechanisms that constitute the overall phenomenon of sediment trans- 
port by ocean waves. After each phase of the problem is well-described 
and understood, the desired relationships will be developed based on all 
available theoretical concepts and empirical evidence. A more detailed 
description of our objective is given in the following section. 


2. FORMULATION OF THE PROBLEM 


It has been observed that loose sediment is moving in considerable 
quantities near and along the ocean floor even in relatively deep water. 
This motion cannot be attributed to the action of ocean currents because 
the flow intensity of these currents is usually very low and the hydro- 
dynamic forces associated with it are not sufficient to overcome the 
forces resisting motion. It is evident, therefore, that the mechanism 
mainly responsible for this motion has its origin in the oscillatory flow 
near the bed which is caused by the surface waves. It is obvious that in 
general there is no net transport associated with this motion since the 
particles oscillate more or less about their mean position. The claim set: 
forth is that the hydrodynamic effect of this flow is to relieve the part- 
icles of all or part of their weight so as to bring them to a state of 
incipient equilibrium. At this state any incidental secondary flow or 
current, no matter how weak, will be able to set the particle in motion. 
The problem now may be defined in a more precise form, An oscillatory 
motion is induced by the surface waves on the boundary layer. It is 
desired 


a. To develop an expression that describes the flow 
field in the boundary layer. 


-b. To determine the dynamic effect of this field on 
the solid particles forming the bed. 


The individual results of these two phases of the problem will be com- 
bined in a logical fashion to obtain a relationship by means of which it 
will be possible to predict the pattern of the motion of solid particles 
near the bed for a given set of wave characteristics and bed composition, 
We proceed first with the study of the flow in the boundary layer. 


3. OSCILLATORY BOUNDARY-LAYER FLOW 
a. Theoretical Considerations 


(i) The Laminar Case. The problem of the boundary-layer flow 
can be treated in the ordinary fashion. According to the fundamental 
principle the flow field outside the boundary layer can be described by 
means of the irrotational theory while the complete equation of motion 
within the layer is simplified with the help of dimensional arguments. 
The boundary conditions of the simplified equation are set so as to 
satisfy the nonslip requirement at the solid boundary and the continuity 
of the velocity components at the outer edge of the boundary layer. Lin 
(1957)* has presented a solution to the more general problem in which the 
flow both inside and outside the boundary layer has a mean steady com- 
ponent as well as an oscillatory one. .In the present case the two flows 
have oscillatory components only. The equations of motion will be 
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inside and outside the boundary layer respectively. Continuity is 
assumed to be satisfied individually 


gu Ov 


so that ae + ay = 0 (3-3) 
and O41 + OM, = 0 (3-4) 
Ox oy 


Equations (3-2) and (3-4) describe a two-dimensional inviscid un- 
steady and incompressible flow. The solution of these two equations under 
a specific set of boundary conditions will define the velocity field in 
the fluid as a function of time and space, provided of course, that the 
effect of viscosity is negligible. In the problem at hand the first 
boundary condition is a sinusoidal progressive wave at the surface (y = 0) 
and the second, zero vertical velocity component at the bottom (y = -d). 
To the first approximation equation (3-2) may be linearized by omitting 
the quadratic terms; it can be written then as 


OU) Sora. MenOy (3-5) 
ot p ax 


* Numbers in parentheses denote date of reference listed on page 29. 


This simplification is justifiable only under the assumption that 
the slope of the surface waves is very small; this is equivalent to 
k = << 1 where k is the wave number and H the amplitude of the surface 
wave. When the water depth is not very large (d<1/k) the solution 
of the simplified equation indicates that the water particles describe 
elliptical orbits around their mean position with component displace- 
ments in the x and y directions given by the expressions: (Lamb, 1932 
sec. 229). 


cosh k (y + d) 


H 
> eran Tel cos (kx-yt) (3-6) 
_H sinh k (y +d) _. 
NOS Crea TT ae ee SESS ed (3-7) 


It is evident that the vertical component of the displacement be- 
comes smaller as the distance from the surface increases and that right 
at the bottom (y = -d) the motion degenerates into a simple harmonic os- 
cillation along the x direction. The corresponding velocity components 
are obtained by differentiation of equations (3-5) and (3-6); so that 


cosh k (y + d) 


Oe ae a = 
uy (x,y,t) = Se ). Srean Tal Sin (kx wt) (3-8) 
A OMe H sinh k (y + d) th ie 
vy Gey tf) = NEE San Ra ee oe (kx - wt) (3-9) 
As y — -d vy (-d) — 0, and 
H : 
u, (x,-d,t) == w cosech kd sin (kx - wt) (3-10) 


il 


or ou, (x,-d,t) = aw sin (kx - wt) = ug sin (kx - yt) 


: H 
where obviously a = 5 cosech kd and ug = ay 


On the basis of dimensional considerations it is reasonable to postu- 
late that the thickness of the boundary layer is very small <6 = iI 
Ww 


compared to 1/k so that for all practical purposes u, (x,t) may be 
assumed constant within the boundary layer and approximately equal to 
u,(x,-d,t); as it is customary we will use the notation U, for the free 
stream velocity at the outer edge, so that 


U,, = Uy (x,-d,t) = aw sin (kx - wt) (3-11) 
with surface waves of large wave length (ka<<l) equation (3-11) becomes 
U,, = aw Sin wt = ug sin wt (3-12) 


Therefore at the edge of the boundary layer equation (3-5) will take 
the form 


oe = i ae (Gaia) 


Subtracting now (3-12a) from the complete boundary layer equation 
(3-1) we obtain 


2 
gu + u gu + V gu = QU co + “ a (3-13) 


ot Ox oy ot oy 
The value of u and v is zero at the wall and equal to U~ and Vy respect- 
ively at distance § from it, where § = fe within the boundary layer, 
therefore to the first appuos ina tn |e = (ee: Moreover it is reasonable 
\ 


to asume that in the same region 


a al HOE Wy) 
oy 5 f ax 
From continuity then V = $k Use 


Let us now examine the order of magnitude of the terms in equation 
(3-13) 
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So we see that the quadratic terms are smaller than the rest by a factor 
ka and since according to our basic assumption ka << 1 these terms can be 
omitted, 


The boundary layer equation (3-1) to a first approximation may be 
written now as 


2 
Ol 5 aula 4, Ga (3-13a) 
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with boundary conditions 


we CO, i) O 


\ (3-14) 


u(o, t) =U = u,sinwt 


The solution of (3-13) satisfying (3-14) is of the form 


W Gy &) us { sinwt - e '’sin (wt-By) } (3-15) 


wo 
2y 


We recognize the second term in the parentheses of (3-15) as the 
solution of the diffusion equation 


where 8 


el ame eso) (3-16) 


with boundary conditions 


u (o,t) Sinwt | } 


iH] 
[= 
fo) 


(3-17) 


‘u (o,t) (0) 


which describes the flow near an oscillating flat plate (Schlichting, 1955). 
The structure of the solution (3-15) suggests that it is possible to deter- 
mine the velocity component in the x direction within the boundary layer 

of an oscillating body of water by means of simple superposition of the 
irrotational component at the bottom and of the solution of the oscillating 
flat plate. This property has significant importance especially in con- 
nection with the experiment. The solution thus obtained is a good ap- 
proximation of the actual case provided that the flow within the boundary 
layer is laminar, 


However, even when the surface of the wall is hydraulically smooth 
one would expect that at a certain value of the Reynolds number defined 
as 


the flow will become unstable. The knowledge of the critical value of 
Reynolds number is very important to our problem, because the dynamic 
effects on the particles under laminar flow conditions are quite different 
than in turbulent flow. It would be very helpful, therefore, to establish 
criteria of flow stability covering a wide range of wave characteristics 
and bed roughness. This is rather an involved case in the class of prob- 
lems of hydrodynamic stability. If we recall the lengthy calculations 
required to solve the relatively simple case of the Blasius profile on a 
smooth flat plate we may conclude without much deliberation that the 
attempt to seek a theoretical solution for our case, of a mean unsteady 
flow and a rough plate, does not offer much hope for success. Besides 


this phase of the problem is beyond the scope of the present study. In 

a practical application it is not as important to know the exact value of 
the critical Reynolds number as to be able to predict with sufficient con- 
fidence that under the existing conditions the flow regime in the boundary 
layer is not laminar and consequently that the theoretical laminar solu- 
tion is no more applicable. This type of information can be obtained by 
experimental methods. The studies of Li (1954) and Manohar (1955) are 

two outstanding sources of such information, The procedure used by these 
two investigators as well as their results will be discussed briefly in 
Appendix A. 


(ii) The Turbulent Case. After we have established the fact 
that the flow in the boundary layer is unstable we proceed with our main 
task which is the description of the velocity distribution within this 
layer. We define the instantaneous velocity components in the two di- 
rections as 


iS mh) 2 me (a1) 


and V=av 


where u is a simple harmonic function of time and space while u' and v' 
are the turbulence components. We claim that the amplitude of u is a 
function of y only and that its frequency is the same as the frequency 
of the wave. A reasonable expression for U will be then of the form (a 
form similar to equation (3-15) of the laminar case) 


@ Gt) =o, {sinwt ~ £1 (y) sin [ wt -£,(9) | } (3-19) 


where fi (y) and f, (y) are functions of y only. These functions can be 
determined by experimental methods, 


This, of course, implies actual measurement of the velocity at points 
very close to the boundary. The effective thickness, however, of the 
boundary layer in the actual case does not exceed a few millimeters and 
even if it were possible to make velocity measurement within such a thin 
layer the experimental flume must have practically prototype dimensions. 
The experimental work can be simplified considerably by making use of the 
principle of superposition mentioned above. Therefore, setting 


au SUS (3-20) 
where Cie U, = us Sinwt 
and u, = uf, oy) sin [ wt-£,0) | (3-21) 


the velocity component u in the boundary layer for given values of a and 
w can be obtained from measurement of f,(y) and f,(y) in a flume in which 
the water surface is maintained at rest while the bottom is oscillating 


with a simple harmonic motion in the x direction prescribed by a and Ww. 
The size of the flume can then be kept within reasonable limits since 
f(y) is a rapidly decaying function of y. 


b. Experimental Work 


Velocity measurements were made in an experimental flume in which 
a horizontal plate near the bottom was oscillated in its own plane with 
a simple harmonic motion with respect to the water above. A detailed 
description of the apparatus and the measuring device has been given 
elsewhere (Kalkanis, 1957). In that previous study experimenting with 
a smooth plate it has been found that the expression describing the 
vertical distribution of the horizontal velocity is of the form 


= € : 
Uy = ugk, (By) sin(wt-a £n koBy) (G22) 
w 4 
where 8 = Ca the values of the exponent c and of the constants as the 


results of these experiments indicated were the following 


ES O65, ty = O82, @ S105, ky = B.888 


The above expression describes the velocity distribution relative to a 
fixed boundary; the distribution relative to the oscillatory boundary 
according to equation (3-20) is given by the expression 


E =u, {sin wt-k,(By)° sin(wt-a én kapy) } (28) 


At the time of the conclusion of the experimental work of this first 
phase of the study the lack of sufficient data led to the adoption of the 
parameter 1/8 as the proper length scale. The data collected in subse- 
quent tests, however, indicated that more suitable characteristic lengths 
were the amplitude of the oscillation a for the case of the hydraulically 
smooth wall and the parameter aQD for the rough wall, where D is the 
roughness diameter (either two or three dimensional). 


The analysis of the new experimental data which can be found in 
Appendix B led to the following expression: 


a al 
= = [2 + £,2(y) - 2f,(y) cos f(y) ii sin(wt+@) (3-24) 
-1 f1(y) sin fa(y) 


where Qcatani oy See eee 
1-f,(y) cos f5(y) 

The form of the functions fj (y) and fo(y) that seemed to best fit the 
experimental data is the following 


752 
(i) Smooth plate f,(y) 0.3 e @ 
1/3 


Z(G) = 2.55 (Gy) 


(ii) Two-dimensional roughness fy) = 28D) 
2/3 
£6) = 0.5 (py) 4 
4h ! Sige 
(iii) Three-dimensional roughness f(y) = O52 agD 
2/3 
fi) = 0.5 (py) i 


We observe that these equations do not satisfy the boundary condition at 
the wall which is not a very serious limitation since it can be easily 
circumvented by assuming the formation of a laminar sublayer as ina 
steady mean flow. 


We want to emphasize, however, that the above expressions are nothing 
more than a convenient representation of actual measurements. It may be 
that the true distribution follows some other law which conceivably will 
be described by an equation derived by a more rigorous analytical process. 
Nevertheless, one will expect that close to the boundary, which is the 
region we are mostly concerned with, the results by the two methods will 
be closely the same since the experimental work was conducted within the 
range of conditions which most likely will materialize in an actual case. 
We claim therefore in concluding this phase of the study that by means of 
the expressions given above it is possible to describe the mean velocity 
distribution near the ocean bed in terms of the surface wave character- 
istics, the depth of the water and the grain size of the uniform material 
forming the bed, 


4. HYDRODYNAMIC EFFECTS OF THE FLOW ON THE BOUNDARY 


The strong interaction between the flow and the loose material form- 
ing the bed of a stream has been studied at great length by numerous 
investigators. The results of these studies led to the development of 
theories describing the phenomenon of sediment transport in rivers. It 
has been reported some time ago, (Einstein, 1948), that a basic similarity 
exists between the motion of sediment in a river and in a large body of 
water in which the main motion is oscillatory. This, of course, does not 
imply that the laws governing the two phenomena are exactly the same, nor 
that the theories describing the motion in one case are directly applic- 
able to the other. It is, however, reasonable to assume that some funda- 
mental concepts used in the derivation of one theory may be applied in 
the derivation of the other. Since we intend to use in our study some of 
the basic principles associated with the motion of sediment by steady 
streams of water we consider it appropriate to give a brief outline of 
the historical development of theories related to this problem. 


a. Theories of Sediment Transport in Rivers 


The older equations which were used to determine rates of sediment 
transport in rivers have been derived mostly by means of empirical methods, 


The basic concept regarding the pattern of motion was that the loose bed 
is sliding in layers under the action of the flow above. The top layer 

of the bed is set into motion by the “tractive force™ or shear, which in 
flows where energy is dissipated mainly to overcome friction, is equal per 
unit area of the bed to the product of the unit weight of the water, of the 
depth of flow and of the energy gradient. When this force becomes larger 
than the force resisting motion per unit area, which is proportional to 
the submerged weight of the particles forming the bed, the latter begin 

to move. The rate of transport determined experimentally was found to be 

a function of the difference between the two forces... No effort has been 
made in developing this theory to explain the actual mechanism of inter- 
action between the solid particles and the flow field. General information 
therefore cannot be deduced from this theory and its application is neces-— 
sarily confined to the narrow range of conditions used in its derivation. 
A major weakness is that it deals only with an average value of the shear 
which is assumed in each case to be constant in time and space. The 
implication is that all particles of a certain size will start moving 
simultaneously over the entire bed, whenever the average shear is larger 
than critical, This, of course, is contrary to the well-established fact 
that motion near the bed takes place in the form of sudden jumps by indi- 
vidual particles alternating with rather long periods of rest. 


Jeffreys (1929) perhaps, was the first investigator to base a theory 
on the stability of the individual solid particle. He suggested the appli- 
cation of the solution from classical hydrodynamics regarding the stability 
of a long circular cylinder of radius r resting on the flat bed of a deep 
stream with its axis perpendicular to the flow. The expression of the 
complex potential is then of the form 


W = mr U coth mfr/z (4-1) 


and the upward thrust transmitted to the solid can be shown to be equal ‘to 


1 2 BO 
= = ied (hm 
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It follows that the condition of motion is given by the inequality 
a. A 2 u 2 a Po-Pr “ (4-3) 
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where Uc is the critical value of the stream velocity. He postulated that 
for three-dimensional elements the values would be slightly larger. The 
numerical values given by Jeffreys for water and sand were Uc = 4.32 cm/sec 
for r = 0.01 cm and U, = 13.6 cm/sec for r = 0.1 cm. This model, of course, 
gives the correct answer as long as there is sufficient justification for 
describing dynamic effects in flows of real fluids by means of the irrota- 
tional theory. Jeffreys answer to this question is affirmative; he claims 
that during the initial stage, when a particle is just dropped on the bed 


of a steady stream, the theory is still applicable because the flow 
around the particle has not yet been modified by viscosity. This may be 
true, but even if we accepted the validity of the argument, still we are 
faced with a phase of our problem which Jeffreys’ model does not seem to 
recognize; this is the erratic fashion by which particles move on the bed, 
According to his theory the flow which is responsible for the forces 
induced on the particle is uniform and steady everywhere on the bed, re- 
sulting ina uniform force field very similar in character to the one 
associated with the tractive force theory described above. Therefore a 
particle that starts moving at some point of the bed will never have a 
chance to come back to rest at some other point on the bed, a mode which 
is inconsistent with the actually observed form of motion. 


The more modern theories go a step further; they concern themselves 
again with the stability of the individual particle, but they recognize | 
the fact that the hydrodynamic forces acting upon it vary rapidly with 
time, a phenomenon strongly associated with turbulence. Moreover, they 
use some elementary concepts concerning the structure of turbulence to 
explain the mechanism of suspension. Accordingly, the particles are being 
moved upwards from lower layers of high concentration to higher ones of 
lower concentration by the vertical velocity fluctuations. Equations have 
been derived based on the concept of momentum transfer by turbulence which 
describe the distribution of concentration of sediment in suspension re- 
lative to a specified value at a reference level. If there were a way of 
predicting this value the rate of sediment moving in suspension would be 
readily obtained by multiplying corresponding values of local mean velocity 
and concentration and integrating over the entire depth. One method to 
predict the concentration at some point in the flow has been proposed by 
Lane and Kalinske (1939). These authors claimed that the bed itself may 
be taken as the reference level. The corresponding concentration can be 
obtained by making use of the statistical properties of turbulence close 
to the bed. This implies direct exchange of particles between the sus- 
pension load and the bed, 


Another method proposed by Einstein (1950) interposes a thin layer 
between the suspension load and the fixed bed. The reasoning behind this 
model is that very close to the bed the scale of turbulence is so small 
that the eddies are of about the same order of magnitude as the particles 
on which they act and, consequently, they are unable to move them away 
from the bed, Within this thin layer adjacent to the bed and called the 
"bed layer" the particles move by sliding and rolling, in a fashion much 
different than the particles in suspension above. The bed layer has been 
assumed to have a thickness of about two grain diameters, an estimate 
based on observation. There is a continuous exchange of grains between 
the bed layer and the bed and between the suspension load and the bed 
layer. The basic concept of the theory of the bed-load function as it 
applies in a river flow is that at equilibrium all these exchanges occur 
at the same rate. The fraction of the total load that is carried within 
the bed layer is called the bed load. The rate of deposition from the 
bed layer to the bed is found to be a function of the bed-load rate, 


while the rate of removal of grains from the bed is a function of the 
local flow intensity. The functional relationship between the "bed-load 
rate" in a stream and the flow intensity constitutes the “bed-load function" 
while the equation expressing this relationship is defined as the "bed- 
load equation". With the help of this equation it is possible to calcu- 
late the bed-load rate for given flow conditions and bed composition. 
Finally, the concentration in the bed layer which can be easily calculated 
from the bed-load rate is assumed to be equal to the concentration of the 
suspension load at the reference level. This helps to find the complete 
solution to the problem which expresses the total rate of transport as the 
summation of the suspension rate and the bed-load rate. The significance 
of the concepts of the bed layer, the bed-load function and the bed-load 
equation is evident. In the absence of substantial evidence to the con- 
trary it is legitimate to assume that similar concepts hold true in the 
case of an oscillatory mean flow. Accepting a priori the existence of a 
bed layer and of a bed-load function, our objective will be to develop 

the bed-load equation associated with this type of flow. In our effort 
along this line the valuable experience from the steady mean flow will be 
used as a guide. The procedure leading to the derivation of the bed-load 
equation will be described in the following section. 


b. The Bed-load Equation 


It is a well-established fact that the distortion of the flow field 
around a solid particle resting on the bed of a stream generates a lift 
force acting on the particle, even if the latter is well-sheltered within 
the sublayer, Naturally the larger the size of the particle relative to 
the thickness of the undisturbed laminar sublayer, the more pronounced is 
the distortion and the greater the intensity of the lift force. This 
force will tend to dislocate the particle and move it away from the solid 
bed. As long, though, as the particle is still in contact with the bed, 
the lift force is acting only vertically upwards and it can be expressed 
as 
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Ci is the coefficient of lift, A, a shape factor and u the instantaneous 
velocity acting at a distance yo from the theoretical bed. In the case 
of a steady mean stream the location of the theoretical bed and the level 
Yo relative to it at which the velocity must be taken have been established 
experimentally by Einstein and El-Samni (1949). They conducted flume ex- 
periments with plastic spherical balls 0.225 feet in diameter placed in a 
steady stream of water. The theoretical bed has been determined as the 
reference level from which distances should be taken so that the measured 
values of the mean velocity would give the best fit to a logarithmic 
distribution. It has been found that the matching of this profile was 
the most satisfactory when the theoretical bed was taken at a distance 
0.20 D below the top of the spherical particle. The distance yo at which 
the velocity should be measured in calculating the 1ift force has been 


obtained simultaneously with the determination of the lift coefficient Ch. 
The method used to this effect consisted of measuring the differential 
pressure in the fluid between two layers; cone at the top of the spheres 
and the other near their base. This pressure difference which is a mea- 
sure of the lift force can be expressed as 
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The analysis of the experimental results revealed that Cy; in the above 
expression had a constant value Cy = 0.178 for a wide range of flow con- 
ditions provided that the average velocity was measured at a distance 
0.35 D from the theoretical bed. Since it is practically impossible to 
determine the value Cy in an unsteady mean flow, it would be reasonable 
to assume that it is about the same as in a steady stream. As a matter 
of fact, as it will be shown shortly, the only assumption that is really 
necessary is that Cy is constant throughout the entire cycle of oscilla- 
tion. Before we proceed any further we wish to summarize the statements 
adopted in this study which have their origin in El-Samni's experiments 
with steady mean flows. 


a. The theoretical bed lies at a distance 0.2 D below the top 
of the grains resting on the fixed bed. 


b. The 1ift force associated with the mean flow can be calculated 
by using the velocity at a distance 0.35 D from the theoretical 
bed, 


c. The 1ift coefficient associated with the mean flow has a constant 
value independent of the Reynolds number. 


If it were not for turbulence it would be rather easy to establish a cri- 
terion of stability similar to Jeffreys’ in the form of the inequality 


L > Ww’ (4-6) 
(W’ is the submerged weight of the particle) 
Gis Ci Pr > A,D > (ps = (ove) gA,D (4-7) 


where by u we mean the amplitude of the velocity calculated from equation 
(3-24) with y = 0.35 D. 


In a turbulent stream, however, all the local flow parameters, and 
consequently the local life as well, vary rapidly with time. An accurate 
description of the temporal variation of the local life by analytical 
methods would be possible only if sufficient information regarding the 
structure of turbulence in the boundary layer were available. Considerable 
_amount of effort is made now by numerous investigators toward developing 
new theories or toward improving existing ones on the subject. The fact 
remains that even if we could wait until the theoretical work had been 
sufficiently advanced to permit practical applications, still some 


approximations would be required since most of the studies are concerned 
with steady mean flows and smooth boundaries. In the meantime and due to 
the lack of more reliable information we are forced to depend again on 
experimental evidence. A possible criterion of stability in a turbulent 
stream could be set as 
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where L' is the turbulent component of the 1ift force. The study of the 
variation of L' with time constituted the last phase of the experimental 
work conducted by El-Samni. By measuring the instantaneous values of 

the lift force exerted by a steady stream of water on the plastic spheres 
mentioned above, he was able to show that L’ behaves like a random vari- 
able having a normal distribution with mean zero and standard deviation 
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flow will be similar although with a different numerical value of oO we 
could proceed as follows: Let p be defined as the probability that a 
particle resting at a certain location in the bed becomes just ready to 
move; this implies that 
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where z of course is a dummy variable. 


The values of No and By, can only be determined experimentally. It 
is evident now why it is not necessary to determine the exact value Ge Cr. 
As long as Cy, is constant its effect will be reflected in the value of Bx. 
We may recall now that u is a periodic function of time which as the ex-— 
periments indicated can be approximated by the right hand side of equation 
(3-24). Consequently y is a periodic function too which means that a more 
appropriate form of equation (4-13) would be 
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The geometric representation of the probability p as expressed by 
equation (4-13a) is shown in the sketch below. 


§(U) 
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It is evident from the form of equation (4-13a) that the turbulent field 
was assumed constant during the entire cycle of the oscillatory mean flow. 
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The number Nj of particles of size D per unit area of the bed surface 
that at any instance become free to move and indeed do move is propor- 
tional to this probability as well as to the total population of similar 
particles per unit bed surface, It is evident that 


NaS 5 (4-14) 
where Aj is a constant. 


The rate of transportation qg which is defined as the rate in dry 
weight at which solid particles move across a section of unit width 
oriented perpendicular to the flow can be expressed as 
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where Wj = Ad? bp Ys, Ys being the dry unit weight of the particles and 
Vi the speed of the particle propagation. It should be remembered that 
the motion of the particles is not continuous, but that it consists of a 
sequence of discrete steps. It will be reasonable, therefore, to express 
this speed of propagation as the ratio of the average distance covered by 
a particle in a step and of the total time required for the completion of 
a full cycle of motion. If we denoted the former by £ and the latter by 
t equation (4-15) after the substitution for Wj becomes 
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The distance 4 as in the case of a steady mean flow may be considered as 
a charactsristic length proportional to the grain diameter, so that g = AL D. 
Moreover by definition t can be written as 


BS By Pe (4-17) 


where ty is the part of the cycle during which the particle is at rest 
(between two consecutive excursions) and t2 is the part during which the 
particle is in motion. Experiments with light-weight coarse material and 
steady mean flows in flumes have demonstrated that in general tz is much 
smaller than ty which led to the conclusion that for all practical pur- 
poses we may assume t = t1. Einstein (1950) has defined tj, which he 
called "the exchange time" as a measure of the time required for the 
replacement of a particle that is just being picked up by the flow at a 
certain spot of the bed by a similar particle that is being brought to 
rest at the same spot. Therefore, one may think of tj as a parameter 
which depends on the properties of the particle and the surrounding fluid 
only and which, consequently, is independent of the local flow conditions. 


The simplest such parameter having the right dimensions is the time re- 

quired for a particle to settle through a distance equal to its diameter 
in the fluid at rest. If the settling velocity is denoted by Vs then we 
can write 


D 
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cy being a constant of proportionality 


A basic characteristic of our model is that the criterion of equili- 
brium is governed mainly by the lift force exerted by the flow on the 
particle, Since the coefficient of this lift force assumes a constant 
value even for very small values of the Reynolds number the settling 
velocity at equilibrium may be written as 
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Combining now (4-18) and (4-19) and making the proper substitutions for 4 
and t equation (4-16) becomes 
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which is the "bed-load" equation; after being rearranged (4-20) becomes 
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This is identical with Einstein's (1950) equation (38) provided that the 
material forming the bed is uniform, The significance of equation (4-21) 
is that it describes the equality between the rate of deposition (left 
hand side term) and the rate of erosion (right hand side). The remarkable 
characteristic of either equation (4-20) or (4-21) is that they indicate 
that the bed-load rate is only indirectly related to the flow intensity 
through the probability p. Eliminating this probability between equations 
(4-13a) and (4-20) we obtain the fundamental relationship between flow 
intensity and bed-load rate. 
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Before we conclude this section, it would be necessary to define a 
little better the average distance of travel £ = AyD. A rather small 
value of p implies that only on a small fraction of the bed surface the 
lift force is strong enough to remove a grain of a given size at any 
instance, Therefore a particle that has been lifted by the flow will 
probably come to rest immediately upon completion of the first step since 
it is very probable that the conditions locally will favor deposition. 

In this case £ = A*;D where A'y is the true constant of proportionality 
between distance of travel in a single step and particle diameter, If, on 


17 


the other hand p is rather large, there is a good chance that the condi- 
tions around the spot on the bed where the particle is coming to rest 
after completing a step of length A';D are not favorable to deposition; 
the result is that the particle will be forced to take an additional step 
of length A*;D. This can be repeated a number of times until the particle 
finally finds a point on the bed where it is permitted to rest. This model 
suggests that the actual distance travelled will be proportional to the 
length of each step (A',LD) times the number of consecutive steps made in 
each realization, This number of steps may be thought of as a discrete 
random variable having a binomial distribution with parameter p. The 
probability that the particle will travel a distance At; D is (1-p) which 
is the probability of failure in the first trial. The probability of 
covering a distance 2A';D is p(1-p), the probability that failure: will 
follow one success, In general the probability of covering a distance 

(n + 1) A';D is p™ (1-p) which is the probability of the first failure 
occurring after n consecutive successes, The expected value of the dis- 
tance covered by the particle in a single realization can be expressed 
then by A;D where 
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We introduce now this expression for A,D in equation (4-20) to obtain 
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Solving (4-23) for p we get 
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We equate finally the right hand sides of equations (4-13a) and (4-24) ta 
obtain the very important relationship between flow intensity and bed-load 
rate in the form of equation (4-27) 
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The practical application of equation (4-27) necessitates the determina- 
tion of the constants Ax and By and of the standard deviation T|, of the 
turbulent lift force. The procedure used toward this end will be de- 
scribed in the following section 


5. DETERMINATION OF Ax, By, and TI, 


Suppose that a rather large set of corresponding values of the para- 
meters ¥ and 6 is available. These values were calculated from equations 
(4-10) and (4-26) respectively. The data used in these calculations may 
have been obtained in the field or in the laboratory but in either case we 
are confident that they are reasonably accurate. Our task then consists 
in selecting a set of values of the parameters to be determined so that ¥ 
and ¢ calculated from equations (4-10) and (4-26) will satisfy as closely 
as possible equation (4-27). The determination of the constants in the 
present study was based on values of ¥ and % calculated from experimental 
data. The description of the equipment and the procedure followed in this 
phase of the study is given in Appendix C. In these experiments three 
different sizes of sand were used and for each size nine runs were made 
with different amplitudes of oscillation and frequencies. The values of 
q were directly obtained from the experiment and they were next introduced 
in equation (4-26) to obtain $ while equation (4-10) was used to calculate 
the corresponding ¥. The latter was then plotted against the former on a 
log-log paper as shown in Figure 8. On the same paper a family of curves 
was plotted also representing graphical solutions of equation (4-27) with 
Noas parameter. The four curves shown on the figure were calculated with 
1/No equal to 1.0, 1.5,.2.0, and 2.5. Both Ax and Bx in these calculations 
were taken equal to unity. The left hand side of equation (4-27) was in- 
tegrated numerically with the simultaneous use of normal error tables. The 
three steps of the procedure used to determine the points of the theoret- 
ical curve were the following. First an arbitrary value of Y was chosen 
which was divided by (cos ¢;)* where gj is the midpoint of one of the 
nine equal intervals dividing the quarter of the cycle, so that ¢, = 5°, 
G2 = 15° etc. The probability p which is approximately equal to the 
average value of Pp; where 
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was next obtained. The rearranged form of equation (4-25) was finally 
used to calculate the corresponding $; so that 


$= — (5-2) 


By examining the shape and relative position of the graphs on Figure 8, 
one may conclude that the translation of the experimental $, Y curve along 
the two axes would bring it to a close approximation with the theoretical 
one constructed on the basis of NA 1.5. The coefficients by which the 


coordinates of the experimental $, Y curve should be multiplied in order 
to achieve this approximation express the sought-for values of the para- 
meters A, and B,. It is evident from the foregoing that 
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We may conclude, therefore, that for a given set of wave characteristics 
and grain size of bed material it is possible to determine the bed-load 
rate qp through equations (4-10), (4-27) and (4-26). 


This is an important result, but by no means the final answer to our 
problem. Because of symmetry the rate in one direction during the first 
half of the cycle will be exactly equal to the rate in the opposite di- 
rection during the second half. The net effect of course will be zero 
steady movement. The question is now as to how can we make a practical 
use of the results obtained so far. The argument advanced is that al- 
though the calculated value of dp does not give a direct measure of the 
amount of sediment that systematically moves in some direction, it can 
be used to determine the number of solid particles per unit area of bed 
surface that at any time are exposed to the transportive effort of any 
incidental flow no matter how weak; this flow is not strong enough to 
dislocate the particles from their state of rest, but once it finds them 
in a state of motion produced by the surface wave it is able to move them 
forward, The rate of sediment transport per unit width of the bed associ- 
ated with such a secondary flow will be 
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where h is the thickness of the layer within which motion occurs, U(y) 
expresses the velocity distribution of the flow and c(y) the concentration 
distribution of the solid particles in motion. The magnitude of h may be 
taken equal to 2D, a customary approximation for flows of this type. 
Another approximation very common with steady mean flows is that c(y) 
within the layer may be considered constant. The velocity distribution 
U(y) of course remains undefined, but for any particular case is presumed 
known. Equation (5-4) then becomes 
2D 
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which will provice the desired answer. to our problem provided of course 
that co is known, This is the phase of the problem where the results ob- 
tained in this investigation find a direct and rather important application, 
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We have shown already that qp can be determined for a given set of con- 
ditions, In the following section a method will be described by means 
of which cg could be calculated from known values of dp - 


6. DETERMINATION OF co 


The rate of flow through a cross section of unit width and height 
h = 2D is 
2D 


Q, 1, UCy) ay (6-1) 
(eo) 


while the rate of sediment transport through the same section at any 
instance is 
2D 
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Co is a measure of the concentration of the sediment which at any time is 
at a state of motion within the bed layer. In this layer the rate of trans- 
port due to the oscillatory flow is by definition equal to dp. Therefore 
in a way similar to (6-2) we can write 
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value of @ which we will denote by 4,. We can write then 
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u is obtained from equation (3-24) and naturally is a function of both 
y and t. Hence the expression for the mean value will be of the form 
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The calculation of 1, from (6-5) is a lengthy and tedious operation which 


has to be performed in each particular case. This renders the method 
proposed: here practically inapplicable. 


A simpler approach could be based on the assumption that Um is pro- 


portional to the amplitude of the velocity at some arbitrary distance 
from the wall within the bed-layer. In other words we postulate that 


Mg =A, [u, | (6-6) 
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with A; a constant of proportionality and Up the velocity from equation 


(3-24) at an arbitrary distance, say y = D, from the theoretical bed, 
A_ now can be calculated from the equation 
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A large number of values could be obtained from (6-7) for a wide range 
of variation of the independent variables, and then averaged out. This 
average value may be considered as a universal constant and be used to 


calculate u, from (6-6) With Up = [ally =n: 


The work could be substantially reduced by the use of a computer; 
otherwise the operation is not much easier than the more accurate one 
mentioned previously. Its main advantage over the latter is that it has 
to be performed only once. Because computer time was not available when 
the study reached this point, a simpler but less accurate method was used 
to calculate u,. This method is adequate at least in establishing an 
order of magnitude, 


The method consisted of a numerical integration of equation (3-19) 
in which the two functions f, (y) and fz (y), as we have seen in section 


3, were of the form 
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One may observe that for given values of the parameters a8 and gD both 
f;(y) and fo(y) are functions of the ratio y/D only. AS a characteristic 
set of these parameters, the mean values from the twenty-seven runs men- 
tioned in a previous section were used, The ratio y/D was made to vary 
in increments of .4, beginning with y/D = .2 and ending with y/D = 1.8. 
The angle wt was varying between 0 and 27 in increments of 7/12 beginning 
with wt = 71/24 and ending with wt = 47 7/24. Typical velocity profiles 
constructed this way and corresponding to different values of the angle 
wt are shown of Figure 10. The value of the constant A5 obtained from 
the integration of these profiles was found to be equal to .618, Equation 
(6-4) can be written now as 


Co = .618 nee (6-8) 


This is the final result of our study. When c, from (6-8) is intro- 
duced in equation (6-2) the desired value of Qp is obtained provided, of 
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course, that the velocity distribution U(y) of the secondary flow is 
known. In actual cases it is hard to predict the character of such 
secondary motions because in general they depend upon local conditions. 
It is, however, possible to deduce reasonable estimates of their behavior 
by statistical methods based on long time records. Since the wave char- 
acteristics themselves are usually evaluated by similar methods it be- 
comes evident that the error in the estimated values of the parameters 
entering our problem as independent variables has a bivariate distribution. 
This, of course, reduces the accuracy of the results. There is at least 
one case, though, in which the steady mean motion U(y) in a body of water 
is induced by the surface wave itself, This particular case will be 
described in the following section, 


7. A PARTICULAR CASE OF THE SECONDARY DRIFT 


This type of a second-order drift which is generated by the surface- 
waves in a direction parallel to the wave propagation has been originally 
studied by Stokes (1851) and more recently by numerous investigators, The 
works of Bagnold (1947) and Longuet-Higgins(1953) wil be singled out 
because in addition to their scientific merits offer a rather convenient 
application to the problem at hand. Associated with this secondary motion 
is a steady mean water particle velocity which usually is called '"mass- 
transport velocity". Stokes' expression for this velocity, which we will 
denote by U, is of the form 
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where y, is measured from the mean free surface negative downwards, The 
only necessary and sufficient condition to be satisfied for the derivation 
of this expression is that the flow is irrotational. Longuet-Higgins 
pointed out that the requirement of small amplitude surface wave, as was 
suggested by Stokes, is not necessary. According to Stokes’ theory, the 
velocity near the bottom is negative (in opposite direction to the wave 
propagation) and for typical values of the product kd it increases with 
distance from the bottom attaining its maximum positive value at the free 
surface. Bagnold (1947) on the other hand has shown experimentally that 
the actual behavior is quite different. He observed a strong forward 
velocity near the bottom and a weaker backward velocity at higher levels. 
This confirmed the belief that Stokes’ theory was not accurate, as it has 
been evidenced from older experiments in which. forward velocities were 
observed both near the bottom and the free surface and backward motion in 
the interior. A more reliable theory which confirms the experimental re- 
sults has been developed by Longuet-Higgins (1953). The two fundamental 
assumptions are that the mean motion is periodic in time and that it can 
be expressed as a perturbation of a state of rest. The theory recognizes 
the existence of three distinct regions, the first near the free surface, 
the second in the interior and the third near the fixed boundary. Associ- 
ated with each of these regions is a particular mode of secondary motion. 
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In the present case our interest is concentrated on the boundary layer 
near the bottom only. Omitting the lengthy procedure of Longuet-Higgins’ 
rigorous derivation it will suffice to present the final results for the 
case of progressive waves in water of uniform depth. The first-order 
motion of the fluid is described by means of the irrotational theory so 
that the horizontal component of the velocity at the boundary may be ex- 
pressed, as we have seen in section 3, as 


e = awSinwt (7-2) 


The "mass-transport velocity" can then be obtained from the second-order 
approximation. 


Ta eae £ (w) (7-3) 
4 4 sinh@kd Pe a 
where f (p) = 5-8 e Mcosu ra) ena (7-4) 


These are Longuet-Higgins"* equations (254) and (253) respectively. 


Consistent with our notation 
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So that the expression for U may be written as 
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where y, measures the distance upwards from the effective bed. It is 
evident therefore that the rate of sediment transport of a certain size 
can be determined by combining the "sediment-transport equation" (6-2) 
and the "mass-transport velocity equation" (7-5). A given set of surface 
wave characteristics, temperature and depth of water and grain size of 
uniform bed material is sufficient for the calculation of this rate. 


8. SUMMARY OF RESULTS 


The independent variables entering the problem of the sediment trans- 
port by wave action are the following: 


The amplitude H, the length ) and the angular velocity w of 
the surface wave. 

The depth d and the temperature T° of the water (and 
consequently its kinematic viscosity y). 

The grain diameter D and the density pg of the granular 
material forming the bed which is assumed to be uniform, 


A method has been described in the preceding sections of this study which 
can be used to predict the rate of sediment transport associated with any 
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given set of values of the independent variables, The procedure to be 
followed consists of the following seven steps. 


a. The necessary condition that has to be satisfied for the method 
to be applicable is that the flow within the boundary layer caused by the 
surface wave is unstable. To test whether or not this condition is satis- 
fied in a particular case the graph in Figure 1 may be used. The roughness 
size ¢ can be taken equal to D while the velocity ug is calculated from 


‘where a= - cosech kd = ; cosech = 


The point with coordinates D and u,/y is plotted in Figure 1 and accord- 
ing to its position relative to the experimental curve for the three- 
dimensional roughness, a prediction can be made about the character of 
the flow in the boundary layer. 


b. Assuming that the test in step a proved that the condition of 
instability is satisfied, we proceed with the determination of the flow 
field in the boundary layer. It has been shown in section 3 that the mean 
unsteady velocity in the layer is given by the expression 


: s 
te = [2 2 9G) = BenG Cos ey GY | sin(wt + @) 
where 9 = feoyat fi(y) sinfg (y) 


1-f,(y) cosf,(y) 


It was also shown that 
_ 133y 


f(y) a 65. af D 


and 
2/3 
ey) = 2) (yy) f 


where eA = OR 
Ny ey 


c. The value of the parameter | y| is next calculated from equation 
(4-10); i.e. 


Og. 7 Pe ep) 
ee cle 
The velocity u is measured at a distance 0.35D from the theoretical bed, 
which means that Bet) 
fi (y) = 5.2 ap 
and 2/3 
Eo(y) = .25 (BD) 
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d, With | Y] known the curve through the experimental points in 
Figure 8 can be used to obtain %. This curve is very closely approximated 
by the graphical representation 0 the equation 


1/5 =Z 2712 
= i in eal oh 
latyy Sais 


2 
cos E 


e. The oscillatory sediment transport rate qp in the bed layer 
which we may call the "oscillatory bed-load rate" is calculated from $ 
through equation (4-26), i.e. 


q fe) = 
epee Disea OO n73/2 
4 B(p.-p,) 


iS) 


f, The concentration Co of the solid particles that at any time 
happen to be in a state of motion within the bed layer (of thickness 2D) 
is q 


= 1 ——— 
Cy .618 aD Gig) 


where up is the value of u from (3-24) at y = D. 


g. The rate of sediment transport per unit width in the direction 
of any incidental secondary flow described by U(y) will be 
2D 


Ory Ree | U(y) dy 
(@) 


In the general case U(y) is an additional independent variable which has 
to be determined by methods similar to the ones used in the determination 
of the surface wave characteristics. In the absence, however, of any such 
flow the only possible steady mean motion in the boundary layer is due 

to the surface wave itself which is called the "second-order drift flow". 
This steady flow within the boundary layer is in the direction of the wave 
propagation, The expression describing the velocity distribution associ- 
ated with this flow is of the form 


= eee nad NIG -By -2By 
U = Oy Sine 5-8 e cosBy+3 e } 


as proposed by Longuet-Higgins. This expression can be substituted for 
U(y) in the above integral to calculate Qg. This is the sought-for value 
of the rate of transport of bed material in the direction of the wave 
propagation. 


9. DISCUSSION OF THE RESULTS 


The procedure outlined in the preceding section leads to a quantita- 
tive answer of our problem. In trying to apply the method to an actual 


26 


case some caution is warranted with regard to its applicability and its 
accuracy. In developing the various relationships it was necessary to 
simplify the physical model by making certain assumptions. Unless these 
assumptions are still valid under the prototype conditions the method may 
need modification. 


The basic assumption which practically governs the entire study was 
that the amplitude of the surface wave was small (kH << 1); this permitted 
the linearization of the equation of motion in the free stream. A second 
assumption, but of lesser importance was that the depth of the water d was 
rather large and uniform (small bottom slope.) 


We wish to emphasize at this point that the proposed method is only 
meant to predict the motion of the bed material within the bed layer, 
Consequently it is not applicable in regions where the mixing is violent 
and where considerable amount of sediment is being carried in suspension 
as for instance near the surf zone and onshore. 


Another point that calls for attention is that the criterion of 
instability of the boundary layer is not well defined. As stated in more 
detail in Appendix A, there is some uncertainty regarding the slope of 
the empirical curve serving as the criterion of instability in the region 

4 -1 : ; : ; : : : 
pa <6x 10 ft  . (Figure 1). Until this uncertainty is removed, it will 
Vv : 
be advisable in this region to use the line with the steeper slope pro- 
posed by Manohar. 


The accuracy of the results, as one would expect, depends entirely 
upon the quality of the approximations made in the course of developing 
the various relationships. We shall examine these approximations of the 
procedure step by step and make suggestions regarding their possible 
improvement. 


a. This step deals only with the applicability of the method and 
was covered already in the discussion above. 


b. The velocity distribution in the boundary layer was determined 
experimentally in a flume under conditions similar but not identical to 
those of the prototype. The error introduced this way cannot be very 
significant and in fact it can be absorbed by the effect of the approxi- 
mation made in the subsequent step. 


c. The value of ¥ was calculated from the velocity profile at a 
distance y = .35D. This distance was so chosen because it has been proven 
correct in dealing with steady mean flows; yet there is no proof that it 
holds true in the case of an oscillatory mean flow also. The accuracy of 
the method may be improved by determining the distance at which ¥ should 
be calculated more precisely. 
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d. The value of § was obtained from the curve in Figure 8. Its 
accuracy depends on how closely the theoretical curve approximates the 
actual one. The exact shape of this curve can be confirmed only after 
a considerable amount of reliable data from the field becomes available. 


e. No approximation was involved in this step. 


f. The determination of the constant in equation (6-8) from the 
experiment with an oscillating bottom is not quite accurate. The error 
due to the value of this constant is not serious however, and can easily 
be reduced as more information from field measurement is obtained. 


g: In order to determine the Sediment, rate an’ the direction of Gehe 
flow the concentration was multiplied by the average discharge in the bed 
layer. This of course is only an approximation of the real mechanism. 

The concentration in the layer actually remains constant in a statistical 
sense only. The individual grains are continuously picked up and dropped 
by the flow. Therefore the solid particleS and the surrounding fluid are 
subjected to random accelerations and decelerations., The effect of the 
resulting inertia forces was assumed to be reflected by the value of the 
empirical constant. However the accuracy of the method may be improved 
considerably by refining the analysis of the force field in the bed layer. 


10. CONCLUSIONS 


Surface waves are responsible for the formation of an oscillatory 
boundary layer near the ocean floor. For a certain range of values of 
the wave parameters the laminar boundary layer becomes unstable. In this 
case solid particles of bed material are brought to a state of incipient 
equilibrium, primarily as the result of the instantaneous hydrodynamic 
lift. The concentration of solid material, in a thin layer adjacent to 
the bed (the bed layer) which at any instance is thus made free to move 
can be determined by the method developed here, The rate of the transport 
of such material in the direction of a secondary flow can be calculated 
by multiplying this concentration with the average discharge due to the 
secondary flow in the bed layer. The accuracy of the method should be 
tested against actual field measurements. 
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APPENDIX A 
STABILITY OF AN OSCILLATORY BOUNDARY LAYER 


It has been shown in section 3 that the wave motion at the outer edge 
of the boundary layer may be approximated by a simple harmonic motion. The 
stability of this layer and its transition into turbulence plays an im- 
portant role in the problem associated with the transportation of sediment. 
In a laminar oscillating stream the predominant force induced by the flow 
on the particle is a tangential drag. When this force becomes larger than 
the force resisting motion the particle will begin to move. A criterion 
of initiation of movement then can be set as 


TAY K (Oz Pe) g tan 8, (A-1) 
where @ is the angle of repose of the particle. The shear stress T at 
any distance from the wall may be determined from equation (3-15). It is 
evident that the maximum value of T occurs at the wall so that 


ar = LBu, (A-2) 


max ™ By ang 
Therefore the criterion of motion in a laminar boundary layer can be set 
as 

K D ee 
vg = [Eiperve) sD tan 61 J Og 


This is Manohar's (1955) equation (19). 


In an unstable boundary layer on the other hand the predominant force, 
as we have seen in section 4, is the hydrodynamic lift which is associated 
both with the mean unsteady flow and the turbulent perturbations. The 
mechanism through which motion is induced is not only different in character 
for the two cases but also in the unstable case it is more intense by several 
orders of magnitude, 


The purpose of Li‘s investigation (Li, Huon 1954) was to determine the 
factors and relationships governing the transition of an oscillatory laminar 
layer over smooth and rough beds. The mathematical model being so complex, 
a theoretical approach was ruled out from the outset. The effort therefore 
was concentrated to obtain empirical results experimentally. The study in 
the laboratory of the boundary layer created by a surface wave would re- 
quire equipment of very large size. It was instead considered expedient 
to investigate the stability of the boundary layer near an oscillating wall 
in a body of water at rest. It is apparent of course that the patterns of 
flow on the prototype and in such a model are by no means identical but 
the reasoning was that the critical values of the governing parameters 


will be very closely approximated by each other. Following is a brief 
description of the apparatus and the procedure used in Li’s study. 


A glass walled flume 12 feet long, 1 foot wide and 3 feet high was 
used. About 4 inches above the bottom of the flume a plate, 78" x 11-1/2" 
was mounted on a series of rollers so that it was free to move in a plane 
parallel to the bottom. Through a series of steel cables and levers the 
plate was connected to a 1/2-HP AC motor which was mounted external of the 
flume. By means of an eccentric arm the rotating motion of the motor was 
converted into a reciprocating one which very nearly approximated simple 
harmonic motion. The motor was equipped with a speed reducer which allowed 
the period of the oscillation to vary from 1 to 150 cpm. By varying the 
eccentricity of the arm the amplitude of the linear motion of the plate 
could be varied from 2 to 18 inches. The amplitude and the frequency of 
oscillation in the experiments were made to vary within ranges so as to 
simulate prototype conditions associated with waves of 0.4 to 60 seconds 
period and 0.5 to 10 feet height in water about 60 feet deep. 


During the testing the water depth was kept at about 2 feet. The 
shortness of the flume and comparatively shallow depth of water tended to 
cause the formation of a standing wave. In order to reduce this wave to 
a minimum a series of vertical baffles together with a set of heavy floats 
were used, 


Three series.of tests were made in this flume. The first, designated 
as the smooth series, was carried out with the surface of the plate waxed, 
The second, designated as the two-dimensional rough series, was carried 
out with the use of half-round wooden strips or steel rods as roughness. 
The third, designated as the three-dimensional rough series, was carried 
out with the use of sand or gravel as a roughness. In all cases each 
individual run was carried out with a uniform roughness. That is to say, 
all roughness elements for any one run were of the same size. 


The procedure was to fix a particular roughness to the plate, choose 
an amplitude of oscillation, and then increase gradually the frequency. 
The type of the flow regime was determined visually by dropping potassium 
permanganate crystals from the water surface to the bottom. As the crystals 
dropped to the bottom they left a trail of dye. Once at the bottom they 
would oscillate with the motion of the platform. Continuing to dissolve the 
crystals would leave a clear back and forth trail as long as the flow was 
stable. The flow was defined as unstable when the trail left by thecrystals 
would break down and mix with layers above. The recorded frequencies and 
amplitudes at the critical limit together with the size of roughness and the 
water viscosity were combined to form a critical Reynolds number. The re- 
sults obtained in these tests as supplemented by Manohar may be summarized 
as follows: 


According to its performance the boundary can be classified as 
hydraulically smooth, as rough and as a transition from smooth to rough. 


The boundary behaves as hydraulically smooth even with two and three- 


dimensional roughness provided that the parameter B, where ¢ is the rough- 
ness diameter is smaller than a certain value; B of course is equal to 


A - The range of B, for the three classifications determined from the 


experiments is the following. 


Two-Dimensional Three-Dimensional 
Roughness Roughness 
Smooth boundary Bo S (87S Be ello 
Transition ONS S. Oa S abot HISSE= Be <0249 
Rough boundary Tela) eS Be sGae) SS Bhs 


(i) For the hydraulically smooth boundary it was found that the 
critical value of the Reynolds number defined as 


was constant and equal to 400. 


(ii) In the transition region the Reynolds number was of the form 


with constant critical value for the two-dimensional and three~dimensional 


roughness; namely 


640 two-dimensional 


Z 
iT] 


104 three-dimensional 


Zz 
I 


(iii) In the rough region the governing parameter was not a Reynolds 
number but the dimensional quantity 


wee On2 


v 


which was found to have the critical values of 


0.2 

UpE 4 

ars FE BoSA x kO two-dimensional 
0.2 

UoE tes 4 , 

= = 1.78 x 10 three-dimensional 


The summary of the results as reported by Manohar are shown in Figure 1% 
The explanation given by these two investigators to the rather singular 
behavior at the rough region was that the transition from laminar to 
turbulent is mainly due to the instability of the flow along the wake 
formed between the individual particles. The characteristic length then 
is the size of the wake which, although a function of the roughness, 
diameter, is not necessarily proportional to it. This explanation, 
although reasonable, is not quite convincing. The writer's reservations 
were confirmed by some observations made during a more recent phase of 
the investigation. These observations were made while measuring the phase 
shift on the vertical in the same experimental flume. The method used 
involved again the use of dye which this time was injected from above 
through a thin brass tube. It was observed that the flow was definitely 


unstable at a value of “° equal to i aclon for a two-dimensional rough- 


v 
ness with ¢ = .0625'. Moreover unstable flow conditions were observed 
while experimenting with three-dimensional roughness of ¢ = .017' when 


the value of Ee) was about equal to Ons These points when plotted in 
Figure 1 demonstrate that the constant value of the critical Reynolds 
number as defined for the transition regime may well be extended to cover 
the rough case too. This implies that in Li*s and Manohar's experiments 
the flow in this region was already unstable before it could be established 
as such from observations. There is a possibility that the crystals and 
the high-density dye solution were confined within the layer of the dead 
water under the theoretical bed. The fluid in this layer oscillated back 
and forth with the plate and as long as the velocities were small it never 
had a chance to spill over the crest of the roughness elements and mix 
with the flow above. So even if the flow were unstable there was no way 
of detecting it. On the other hand, one may argue that the more recent 
observations do not describe the behavior of the real model either, be- 
cause the observed premature transition might have been triggered by the 
tube itself. We want to make clear, however, that the tip of the tube 
was maintained at a level several millimeters above the bed at a distance 
where the amplitude of the fluid velocity was negligible while the first 
indication of instability was observed near the bed. In concluding this 
Appendix we would like to point out that this phase of the problem needs 
further investigation so that more reliable information will become 
available. 


* Figure number refers to figures following main text. 


APPENDIX B 
VELOCITY DISTRIBUTION IN OSCILLATORY BOUNDARY LAYER 


This section contains the description of the experimental method 
and the results obtained in the phase of the study which had as its 
objective the derivation of an expression describing the velocity dis- 
tribution in an oscillatory boundary layer. The linearized equation of 
motion from classical hydrodynamics may be used to describe the flow field 
in the interior, Within the boundary layer as long as the flow is laminar 
the velocity field can be described to the first approximation by an ex- 
pression of the form 


u = Uo {sinut = oy sin Cwt-By) } (B-1) 


where ug is defined as the amplitude of the velocity at the wall from 
irrotational theory. 


When the flow in the boundary layer becomes unstable the only way 
the field can be described is by empirical methods. The experimental work 
involved may be extremely difficult; however, it becomes much simpler if 
one postulated that the unsteady mean flow in this case can be described 
by an expression of the form 


u = ug 1Sinwt-f,(y) sin | wt-f5(y) (B-2) 
fof 1 . 2 
in analogy with the stable solution. 


Equation (B-2) may be written as 


u u UD | 
= S|) 2) sanat =|= | sin | wt-£0) | (B-3) 
Uo fe) Uo 
font a { 
where ball | =a! and 2 | = f1(y) 
Uo Uy 


The purpose therefore of the tests was to determine the functions f(y) 
and fo(y). The structure of equation (B-2) suggested the possibility of 
determining these functions by the convenient method of measuring the 
velocity amplitude and the phase shift in a flume where the bed oscillates 
with an harmonic motion of the form uy = upsinwt relative to the water 

at rest. 


The flume and the rest of the equipment used in these tests has been 
described in detail elsewhere, Li (1954), Manohar (1955), Kalkanis (1957). 
In Appendix A we gave a brief description of the main apparatus and in 
this section we will give the necessary information regarding the opera- 
tion of the additional equipment used in the tests. 


The two functions f,(y) and f,(y) can be determined independently 
and indeed in some cases they were so. However, the device used in the 
measurement of the velocity amplitude at some distance from the oscillat- 
ing wall was simultaneously measuring the phase angle of the velocity at 
this level relative to the wall. In other words the same record could be 
used to determine both fj(y) and fo(y). This instrument is essentially a 
modified pitot tube shown diagrammatically in Figure 11. The basic prin- 
ciple of its operation is that the differential pressure at the two tips 
causes the diaphragm to deflect; in doing so it modulates a high-frequency 
signal by altering the capacitance of a Rutishauser pressure pickup. The 
modulated signal is rectified by a discriminator into a voltage which is 
proportional to the differential pressure at the tips of the instrument 
and consequently a function of the incidental velocity. Therefore the 
excursions of the needle of an oscillograph which is activated by this 
voltage will give a measure of the local velocity. The instrument was 
calibrated dynamically by oscillating it with a prescribed simple har- 
monic motion in water at rest and correlating the amplitude of the needle 
excursions with the amplitude of the velocity of oscillation. 


The tests were made on a smooth plate of a high-gloss finish as well 
as on rough plates with two and three-dimensional roughness fixed on them. 
The vertical distances recorded during the tests were measured from the 
crest of the roughness elements. These values were subsequently corrected 
to account for the distance between this level and the theoretical bed. 
The correction in all cases was equal to 0.2D. The values of the function 
fi) resulting from the analysis of the experimental records are listed 
on Tables I, II, and III. The graphical representation of the results is 
shown in Figures 2, 3 and 4, The exponential dependence of the velocity 
on the distance seems to describe the data better than any other simple 
relationship. A small number of points only, obtained in some character- 
istic runs, were plotted in order to make the graphs readable. However, 
for the selection of the most representative relationship the complete 
data was used, 


The experiments with smooth wall indicated that the proper character- 
istic length for the normalization of the argument of f(y) should be the 
amplitude of the oscillation. In_a previous report (Kalkanis 1957), it 
was suggested that the parameter = should be used for this purpose in 


analogy to the laminar case. The original work on which this suggestion 
was based consisted of three runs only reported here as runs 111, 112, 
113, a very inadequate number. As part of the present investigation, 
which has as its main objective the study of the rough case, four more 
runs were made with a smooth wall (114, 115, 116, 117) in order to test 
the reliability of the instrument. The second set of measurements was 
made more than a year after the first and considering the fact that the 
equipment in the meantime underwent certain modifications, the agreement 
of the results between the two sets was remarkable. In any case the 
adoption of the amplitude as the characteristic length was strongly sup- 
ported by the experiments with a rough wall performed more recently. 


By adopting a similar purely empirical approach the analysis of the data 
obtained in these tests revealed that the more suitable length scale to 

be used with f,(y) was the parameter a8D. The forms that most closely 

approximated the experimental data were 


= 2 sz 108 
fi(y) = eee for two-dimensional roughness 
and _ 133y 
181 GY) = soe ep for three-dimensional roughness. 


It is understood that we do not claim that these expressions describe the 
real physical model; on the other hand we believe that the flow field 
described by these equations cannot be very much different than the 
actual one since a considerable amount of data was used in their deriva- 
tion. The implication is that the calculated values by means of these 
functions will be as close to reality as the form of equation (3-15) 
permits. 


In the determination of the function fo(y) the characteristic length 
for all conditions of roughness seemed to be the parameter 4 of the 
laminar case. The phase angle seemed to increase as a power function 


with distance. The empirical relationships obtained as the results of 
these tests are the following: 


5S) (ayn Smooth wall 


f5(y) 


and 


55) (any Two or three-dimensional 
rough wall. 


f,(y) 


The phase angle was measured by two different methods; in the first it 
was obtained directly from the velocity measurement record. The time at 
which the flow was changing direction was registered on the time scale 
of the velocity record while the time of change in direction of the plate 
motion was recorded by another needle of the oscillograph. The time 
interval between these two changes of direction after having been aver- 
aged out for a number of periods gave the desired phase angle. In the 
second method the record of the change of flow direction was obtained 
visually; through a thin brass tube a dye streak was introduced into the 
flow; the change of direction of the streak at any level was observed 
and it was recorded by means of a push button arrangement; the motion 
of the plate was recorded as above. This method is more accurate be- 
cause the equipment it employs is very simple; but on the other hand it 
is subject more to personal bias. The two methods were checked against 
each other in tests with smooth wall and the good agreement of the re- 
sults was a proof that either one may be used with equal confidence. The 
first method has been used to measure the phase shift in the smooth case 


and with a wall having a two-dimensional roughness, while with sand as 
roughness the phase shift was measured by the dye method. 


With both f1(y) and fa(y) determined equation (B-3) may now be 
written as 


5 1/2 
I+ tqGy) n= 2tialGy) cos £2(y) | sin(wt-6) 


(B-5) 
As paar [ fi(y) sin fa(y) 
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APPENDIX C 


DETERMINATION OF A,, B, AND 1], 


The method used for the determination of A,, By, and Ne has been 
outlined in section 5 of the text. What we intend to do here is to 
describe the experimental procedure used in connection with this method. 


The same flume which has already been described in previous sec- 
tions was used, but with some slight modifications. The method called 
for actual measurement of the quantity qg. By definition qp is the rate 
at which sediment near the bed crosses a section of unit width positioned 
perpendicular to the wave propagation. 


Our experinental procedure was based on the proposition that the flow 
field and consequently the magnitude of dp when the bed is fixed and the 
boundary-layer flow is caused by a surface wave are the same as when the 
bed is oscillating with a simple harmonic motion while the water surface 
is at rest. The validity of this argument will be discussed later but 
for the time being let us assume that it is correct. The oscillating 
plate had to be modified so that measurement of q,g could be made possible. 
The plan view and longitudinal section of the modified plate is shown in 
the sketch below 


The procedure is self-explanatory. The space between the trays was 
filled with loose sand of a given size. The plate then was set into 
motion and after the completion of a number of periods it was stopped, 
and the amount of sand collected in the middle tray was removed and 
measured, A wire screen was placed on top of the trays so as to be 
flush with the rest of the sandy bed to achieve uniform roughness 


conditions. Three sizes of sand were used each with three different 
amplitudes of oscillation. Thus 27 runs were made altogether. The 
size range used was 


psduel cs) 1D) < #10 mesh 
#10 < OD < # 8 mesh 
#38 < OD << # 6 mesh 


The average value of D for each of these ranges were 5.51 x On feet, 


ToS Se one feet and 9,25 x Oe feet, respectively. As a matter of 
fact, these were exactly the same sizes used in the measurement of the 
velocity and phase shift distribution. The size of the sample in each 
run was measured volumetrically. No distinction was made as to whether 
the material in the tray came over the rizht or left edge. Each run was 
repeated a number of times and the measured quantities were averaged out. 
The bulk volume of the sample was converted into dry weight by multi- 
plication through a coefficient determined experimentally. 


Since the width of the flume was equal to 1 foot, the dry weight 
of the sample when divided by the period and the number of oscillations 
completed during the testing time gave the value of dg. Next equation 


(4-25) was used to determine 6. This equation is of the form 


ate q3 Pe 33/2 
=) a joe SENT) 
Vs &(ps-Pf) 


qB 32. 


p 
For given values of ae the ratio TF is proportional to D Herein 
Ps bee op 3/2 
OE was equal to 2.63 which means that = = 1190 D . the graphical 


representation of this equation shown in Figure 9 was used to calculate 


$ from measured values of qp. 


On the other hand ¥Y was calculated from equation (4-10) which is 


Ps—P D 
Wife) aeons + 
Pr a 


In fact, the amplitude (icilg was used instead of a because what we were 
really interested in was the value of |v]. A pair of values of the 
parameters $ and ¥ were thus obtained for each run, These values were 
next plotted against each other as shown in Figure 8. Although there 

is considerable scatter, it seems that the experimental points can be 
reasonably represented by the empirical curve drawn througn them. We 

may claim therefore that this curve expresses the functional relationship 


between 6 and ¥Y. On the same figure a family of theoretical curves was 
drawn as it has been explained in section 5 of the text. These curves 
are the graphical representation of equation (4-27) in which the values 
Ax and Bx were equal to unity. Each curve corresponds to a particular 
value of ‘Ilo ° One may observe that the theoretical curve with Ss 1.5 


is very similar to the experimental one but offset both in the Heretical 
and the horizontal. A parallel translation of the latter by a factor 30 
along the horizontal and by a factor 4 along the vertical makes it prac- 
tically to coincide with the former. This means that the values of 6 
and Y, satisfying equation (4-27) with Ax = 30, Bx = 4 and No SiS 


under any set of experimental conditions will be very close to the values 
of the same parameters calculated directly from measured quantities. In 
conclusion we set forth the claim that a theoretical curve constructed 
from equation (4-27) with the constants as determined here could be used 
to describe the relationship between the two functions $ and ¥, 


Several reviewers in the past and more recently P. Lhermitte (1961) 
are rather critical of Li's and Manohar's experiments. More specifically 
they express some reservations about the applicability of the results in 
an actual case on the ground that in the experimental flume an inertia 
force is induced on the particle which in the prototype is absent. Their 
contention is that this force is significant and consequently cannot be 
ignored. Since the same flume was used in the present study one may 
anticipate similar criticism especially in connection with the determina- 
tion of qp from direct measurement. 


We believe, however, that the criticism is not fully justified be- 
cause this inertia force under the average experimental conditions is 
indeed small compared to other forces acting simultaneously on the 
particle. The maximum value of the angular velocity in the set of runs 
with a = 1.25 feet was w = 1.86 rad/sec. corresponding to a maximum tan- 
gential acceleration of 


aw? = 1.25 x 1.862 = 4,32 ft/sec? 


which is much smaller than g. The tangential force could have some 
effect in setting the particle into motion if it were in phase with the 
lift in which case the conbined effect could not be ignored. Since the 
two forces are 90 degrees out of phase the instant one reaches its maxi- 
mum value the other is practically equal to zero. Therefore it is 
justifiable to base the condition of equilibrium on the balance of the 
vertical forces whose absolute value is relatively large and ignore the 
effect of a much smaller horizontal component which is fully out of 
phase, 
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